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THE  HANNAY-BERRY  PHASE  OF  THE  VIBRATING  RING 

GYROSCOPE  * 

SEAN  B.  ANDERSSON  t  AND  P.S.  KRISHNAPRASAD  * 


Abstract. 

In  an  analysis  published  in  1890  G.H.  Bryan  investigated  the  retrograde  precession  of  the  nodal 
points  in  a  vibrating,  rotating  shell  and  wrote  down  a  formula  relating  the  rate  of  precession  to  the 
rate  of  rotation.  This  effect  has  been  utilized  in  the  design  of  various  vibratory  gyroscopes  includ¬ 
ing  modern  MEMS-based  devices.  Existing  analyses  model  these  systems  with  a  pair  of  harmonic 
oscillators  coupled  through  the  Coriolis  force  (the  normal  mode  method).  In  this  work  we  utilize 
the  theory  of  moving  systems  developed  by  Marsden,  Montgomery,  and  Ratiu  to  show  that  the 
nodal  precession  can  be  understood  as  a  geometric  phase  with  respect  to  the  Cartan-Hannay-Berry 
connection.  This  approach  allows  us  to  explicitly  characterize  the  simplifications  of  the  linearizing 
assumptions  common  to  previous  analyses.  Our  results  match  those  of  Bryan  for  small  amplitude 
vibrations  of  the  ring.  We  use  the  inherently  nonlinear  nature  of  the  moving  systems  approach  to 
calculate  a  (small)  correction  to  the  rate  of  precession  of  the  nodes. 


Key  words,  geometric  phase,  holonomy,  fiber  bundle,  vibratory  gyroscope,  MEMS,  rotation 
sensing 

AMS  subject  classifications.  53C29,  70G45,  70H05,  70K65,74F99,  74K10,  81Q70,  93C95 

1.  Introduction  .  In  1890  G.H.  Bryan  published  a  paper  on  the  nature  of  the 
beats  generated  when  a  vibrating  shell  is  rotated  about  its  central  axis  [8].  The 
phenomenon  he  describes  is  quite  easy  to  observe;  simply  take  a  wine  glass,  strike  it 
to  produce  a  clear  tone,  and  then  rotate  it  about  its  stem  to  produce  audible  beats. 
Bryan  noticed  that  these  beats  are  the  result  of  a  precession  of  the  nodal  points  with 
respect  to  the  shell  itself  and  provided  the  following  reasoning.  Consider  a  ring  or 
cylinder  rotating  about  its  central  axis  and  vibrating  with  nodes  at  B,D,F,  and  H  as 
indicated  in  the  left-side  image  of  Figure  1.1.  The  material  points  at  A  and  E  are 
moving  towards  the  center  O.  This  increases  their  actual  angular  velocity  above  that 
of  the  imposed  rotation  and  gives  them  a  relative  angular  acceleration  in  the  direction 
of  rotation  as  represented  by  the  arrows  at  A  and  E  in  the  right  side  image  of  Figure 
1.1.  Similarly,  the  material  points  at  C  and  G  are  moving  outwards  and  thus  their 
angular  velocity  is  reduced.  Those  at  B  and  F  are  moving  with  greater  total  angular 
velocity  than  the  rest  and  thus  experience  a  relative  outwards  acceleration  due  to  a 
greater  centrifugal  force.  Finally,  the  material  points  at  D  and  H  are  moving  with  the 
least  angular  velocity  and  thus  experience  a  relative  acceleration  inwards.  Comparing 
the  arrows  in  the  two  images  of  Figure  1.1  reveals  that  the  effect  of  these  relative 
accelerations  is  to  cause  retrogade  motion  of  the  nodes  relative  to  the  ring.  Using 
classical  variational  techniques  Bryan  derived  a  linearized  partial  differential  equation 
describing  the  behavior  of  the  system,  found  a  formula  for  the  rate  of  precession  and 
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Fig.  1.1.  Nodal  precession  in  cylinder  or  ring  (Figure  from  [8]) 


discovered  that  this  rate  is  proportional  to  the  rate  of  platform  rotation. 

Due  to  the  immense  number  of  potential  applications,  research  on  gyroscopes  has 
been  active  for  many  years.  Devices  have  been  proposed,  analyzed,  and  produced 
using  a  variety  of  materials  and  techniques.  Because  they  lack  rapidly  spinning  parts, 
have  low  power  requirements,  and  are  inherently  scalable,  vibratory  gyroscopes  have 
become  particularly  popular  [21].  One  of  the  most  successful  initial  designs  was 
Delco’s  Hemispherical  Resonator  Gyroscope  (HRG)  [15]  due  to  Loper  and  Lynch, 
which  was  able  to  achieve  performance  levels  equal  to  the  best  ring  laser  gyroscopes. 
This  design,  shown  diagrammatically  in  Figure  1.2,  consists  of  a  quartz  hemispherical 
resonator  supported  on  a  central  stem  and  contained  inside  an  evacuated  housing. 
As  predicted  by  Bryan,  the  nodal  points  of  the  vibration  in  the  hemisphere  precess 


Fig.  1.2.  Delco  Hemispherical  Resonator  Gyroscope  (Image  from  [If]) 


with  respect  to  the  shell  as  the  device  is  rotated.  The  HRG  is  driven  into  elliptical 
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vibration  as  shown  in  Figure  1.3  and  the  resulting  precession  rate  is  about  0.3  of 
that  of  the  imposed  rotation.  The  device  operates  over  a  wide  temperature  range,  has 
high  operating  acceleration  ranges,  low  accleration  sensitivies,  and  negligible  magnetic 
sensitivity  [16].  Similar  ideas  were  used  in  the  design  of  a  vibrating  disc  gyroscope  by 


Fig.  1.3.  Precession  of  Delco  HRG  (Image  from  [If]  following  [15]) 

Burdess  and  Wren  [10]. 

With  the  explosion  of  MEMS  technology  constant  innovations  are  resulting  in 
smaller,  cheaper,  and  more  accurate  devices.  Existing  MEMS-based  devices  include 
tuning-fork  [5]  and  vibrating-ring  designs  [4,  22]  such  as  the  one  shown  in  Figure 
1.4  (provided  by  Douglas  Sparks  of  Delco  Automotive  Systems).  Additional  designs 
proposed  include  a  vibrating  cylinder  [9,  28]  and  a  surface  acoustic  wave  generator 

[27]. 


Fig.  1.4.  MEMS  Gyroscope  (Image  courtesy  of  Douglas  Sparks  of  Delco  Automotive  Systems) 


These  gyroscopes  all  take  advantage  of  the  same  physical  effect,  the  Coriolis  force 
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arising  from  the  non-inertial  character  of  the  rotating  frame  of  the  system  [14].  Mod¬ 
ern  analyses  are  linear  in  nature  and  view  the  Coriolis  force  as  providing  a  coupling 
between  two  vibratory  modes  of  the  structure  [25] .  It  is  desirable,  however,  to  have  a 
method  which,  at  least  in  principle,  can  be  extended  to  a  nonlinear  theory  and  which 
provides  a  unified  setting  for  understanding  a  variety  of  systems  in  which  the  Cori¬ 
olis  force  plays  a  role.  Motivated  by  these  considerations  we  are  led  to  an  approach 
developed  by  Marsden,  Montgomery,  and  Ratiu  based  on  modern  developments  in  ge¬ 
ometric  mechanics.  This  method  is  known  as  the  moving  systems  approach  [18].  The 
technique,  descending  from  the  classical  work  of  E.  Cartan  [11],  describes  the  effect  of 
imposed  motion  on  a  system  as  a  geometric  phase  with  respect  to  a  particular  Ehres- 
mann  connection  called  the  Cartan-Hannay-Berry  connection.  This  geometric  phase 
is  termed  the  Hannay-Berry  phase.  In  previous  work  we  have  used  this  approach 
to  determine  that  the  Hannay-Berry  phase  in  a  rotating,  equal-sided,  spring-jointed, 
four-bar  mechanism  is  zero  [3].  (See  also  [2].) 

Interest  in  the  effect  of  geometric  phases  in  physical  systems  was  spurred  by  work 
of  Michael  Berry  in  the  early  1980’s  [6,  7].  Berry’s  initial  result,  surprising  in  its 
simplicity,  has  sparked  a  great  deal  of  research  that  continues  to  this  day  (see,  e.g. 
the  compendium  volume  of  Wilczek  and  Shapere  [24]). 

While  existing  techniques  have  proved  to  be  effective,  as  evidenced  by  the  various 
devices  constructed  from  those  principles,  it  is  to  be  expected  that  a  deeper  under¬ 
standing  will  emerge  by  appealing  to  a  nonlinear,  geometric  approach  directed  at  more 
accurate  constitutive  models.  The  fact  that  Bryan’s  result  is  geometrical  in  nature 
is  intriguing  and  suggests  that  these  modern  tools  may  prove  useful.  In  this  paper 
we  begin  this  study  by  applying  the  moving  systems  approach  to  the  vibrating  ring 
and  show  that  Bryan’s  result  can  be  interpreted  as  the  geometric  phase.  The  main 
contribution  of  this  work  is  to  provide,  through  the  use  of  geometric  techniques,  a 
clear  understanding  of  the  role  of  the  linearization  assumptions  common  to  previous 
analyses  and  to  make  explicit  the  simplifications  that  result.  Furthermore  we  calcu¬ 
late  the  effect  of  the  neglected  terms  by  deriving  a  formula  for  a  correction  term  to 
the  rate  of  precession  of  the  nodes  of  vibration.  Finally,  this  work  illustrates,  through 
the  reasonably  complex  example  of  the  ring  gyroscope,  how  the  methods  of  modern 
geometric  mechanics  are  useful  in  understanding  physical  effects  in  systems  which  are 
useful  for  sensing  and  as  a  result  hopefully  motivates  new  applications  of  the  ideas. 

The  remainder  of  the  paper  is  organized  as  follows.  In  the  next  section  we  provide 
a  brief  background  the  proper  mathematical  setting  for  geometric  phases,  namely 
that  of  connections  on  fiber  bundles.  In  Section  3  we  summarize  the  moving  systems 
approach  of  Marsden,  Montgomery,  and  Ratiu.  In  Section  4  we  apply  the  technique 
to  the  vibrating  ring  and  derive  the  Hannay-Berry  phase  which  is  then  compared  in 
Section  5  to  the  results  of  Bryan.  The  imposed  rotation  on  the  ring  gives  rise  to 
nonlinear  terms  in  the  equations  of  motion  even  when  a  linear  approximation  of  the 
nominal  (non-moving)  ring  is  used  and  in  Section  6  we  calculate  a  correction  term 
incorporating  these  effects.  We  conclude  in  Section  7  with  a  few  comments  on  the 
results  and  on  future  work. 

2.  Background  .  The  natural  mathematical  framework  from  which  to  approach 
geometric  phases  is  that  of  connections  on  fiber  bundles.  In  this  section  we  present 
a  review  of  the  necessary  background  on  fiber  bundles,  connections,  and  holonomy. 
Additional  references  for  this  material  include  [12,  19]  and  references  therein. 

2.1.  Fiber  bundles  and  Ehresmann  connections.  A  fiber  bundle  is  de¬ 
fined  as  follows.  Let  P,  F,  and  B  be  manifolds  referred  to  as  the  total  space  (or 
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bundle  space),  the  fiber,  and  the  base  space  respectively.  Let  tt  :  P  — >  B  be  a 
surjective  submersion.  We  require  that  P  be  locally  a  product  space,  that  is,  for  every 
b  £  B  there  is  a  neighborhood  U  of  b  such  that  7r-1([/)  is  diffeomorphic  to  F  x  U. 
The  fiber  over  b,  7r-1(&),  is  a  diffeomorphic  copy  of  the  fiber  F  for  every  b  £  B.  See 
Figure  2.1.  The  bundle  is  denoted  by  the  triple  (P,F,B)  or  by  the  projection  map 


Fig.  2.1.  Fiber  bundle,  fiber  over  b,  and  splitting  of  the  tangent  space 


7r  :  P  — ►  B.  If  the  bundle  is  globally  a  product  bundle,  P  =  F  x  B,  then  it  is  called  a 
trivial  fiber  bundle. 

Given  p  £  P,  there  is  a  natural  subspace  of  TpP  (the  tangent  space  to  P  at  p) 

called  the  vertical  space  at  p,  denoted  by  Vp  and  defined  by  Vp  =  kerTp7r.  Here  Tptt  is 
the  linearization  of  the  projection  map  evaluated  at  p.  The  union  of  these  subspaces 
over  all  p  is  called  the  vertical  subbundle  V,  i.e.  V  =  U pepVp. 

Definition  2.1.  An  Ehresmann  connection  A  on  P  is  a  vertical-valued  one- 
form  on  P  satisfying: 

1.  Ap  :  TpP  — >  Vp  is  a  linear  map. 

2.  Ap  is  a  vertical  projection.  That  is,  Ap(v)  =  v  \/v  £  Vp. 

■ 

The  connection  defines  a  horizontal  space  F[p  =  kerHp  at  each  point  p  £  P.  The 
conditions  in  the  definition  imply  TpP  =  Vp  ©  F[p  and  thus  the  connection  gives  a 
splitting  of  the  tangent  space  at  each  point  p  into  a  vertical  and  a  horizontal  part 
(see  again  Figure  2.1).  The  union  of  these  subspaces  over  all  p  is  called  the  horizontal 
subbundle,  i.e.  H  =  U p^pHp. 

Given  an  Ehresmann  connection  A,  a  point  p  £  P,  and  a  tangent  vector  w  £ 
T7T(p)B,  define  the  horizontal  lift  of  w  to  TpP  as  the  unique  tangent  vector  in  Hp  that 
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projects  to  w  under  Tptt.  We  call  this  lift  horp.  The  lift  of  w  can  be  found  by 

(2.1)  hoTpw  =  w  —  Ap(w) 

where  w  £  TpP  is  an  arbitrary  vector  satisfying  Tptt(w)  =  w. 

Lemma  2.2.  The  map  horp  is  well-defined. 

Proof.  Let  ui,  u2  £  TpP  be  tangent  vectors  such  that 

Tp7t(ui)  =  Tpir(u2)  =  U. 


Notice  that 


Tp 7r  (ui  —  u2)  =  Tptt  (hi)  —  Tpn  (u2)  =  u  —  u  =  0 


and  thus  u±  —  u2  is  a  vertical  vector.  Let 


Then 


hor*  [u]  =  iii  —  Apui 


and  horp[u]  =  u2  —  Apu2. 


horp  [«]  —  horp  [it] 


(Ui  -  ApUl)  -  (Sis  -  Apu2) 
(Ui  -  u2)  -  Ap  («1  -  Ui) 
(ui  -  Ui)  -  (mi  -  Ui)  =  0. 


□ 

It  should  be  noted  that  while  here  we  have  defined  the  Ehresmann  connection  as 
a  vertical-valued  one-form  and  derived  the  horizontal  space  and  horizontal  lift,  one 
can  also  begin  with  the  definition  of  the  horizontal  space  or  the  horizontal  lift  and 
define  the  other  two  objects.  See  [18]  for  details. 

2.2.  Parallel  transport  and  holonomy.  Given  an  Ehresmann  connection  A 
on  a  fiber  bundle,  define  parallel  transport  with  respect  to  A  along  a  curve  lifted  from 
the  base  space  in  the  following  way.  Let  6(f),  t  £  [0,1],  be  a  piecewise  differentiable 
curve  in  B.  The  horizontal  lift  of  6(f)  with  respect  to  A  is  the  curve  p(t)  in  P  such 
that  7T (p(t))  =  6(f)  and  that  the  tangent  vector  is  horizontal  for  each  f  G  [0, 1], 
We  have  the  following  proposition  from  [18]. 

Proposition  2.3.  [18]  Given  a  curve  6(f),  t  £  [0,1],  in  B  and  a  point  po  £ 
7r_1(6(0)),  there  exists  a  unique  locally  defined  horizontal  liftp(t)  ofb(t)  to  P  satisfying 
p(0)  =  po  if  P  is  a  locally  trivial  fiber  bundle. 

Proof.  Since  P  is  locally  a  trivial  bundle  we  can  write  locally  p  =  (6,  /)  and 
p  =  (it,  v)  for  6  £  B,  f  £  F,  u  £  T^B,  and  v  £  TfF.  We  can  then  write  the  connection 
one  form  as 


A(b,  f)(u,  v)  =  (0,  v)  +  A{b ,  f)(u,  0) 

=  (0,  v  T  A(6,  f)u). 

(u,v)  is  horizontal  if  and  only  if  A(b,  f)(u,v)  =  0  and  thus  for  a  horizontal  tangent 
vector  v  =  —  A(6,  f)u.  If  6(f)  is  a  path  in  B  denote  p(t)  =  (6(f),  /(f))  where  /(f)  is  the 
solution  to  the  ordinary  differential  equation 

(2-2)  ^(f)  = -A(6(f),/(f))6(f),  /(0)  = /q  where  po  =  (60,/o). 
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Then  by  local  existence  and  uniqueness  for  ordinary  differential  equations  this  defines 
f(t)  and  thus  p(t)  for  small  t.  If  p(t)  can  be  extended  for  all  t  G  [0, 1]  the  connection 
is  called  complete.  □ 

Given  any  curve  b(t)  in  B,  t  G  [0, 1],  and  an  Ehresmann  connection  A,  the  parallel 
transport  operator  v,  is  defined  as 

(2.3)  Tb  :  tt_1(6(0))  ->  7T—1  (6(1)) ,  Tb(p(0))  =  p(  1) 

where  p(0)  G  7r_16(0)  and  p(t)  is  the  horizontal  lift  of  6(f)  with  respect  to  A  starting 
at  p{ 0). 

By  the  uniqueness  of  the  horizontal  lift,  r&  is  a  bijection  from  7t_1(6(0))  to 
7r_1(&(l))  and  by  the  smooth  dependence  of  solutions  of  ordinary  differential  equations 
on  initial  conditions  it  is  a  diffeomorphism. 

Let  bo  be  an  arbitrary  point  of  B  and  let  Cbo  be  the  set  of  all  closed  curves  at  b0, 
that  is  all  b(t)  such  that  6(0)  =  6(1)  =  bo-  The  diffeomorphism  of  7r^1(6o)  onto  itself 
given  by  parallel  transport  along  6(f)  is  called  the  holonomy  of  the  path  6(f).  Let 
<1>60  be  the  collection  of  all  parallel  transport  operators  over  Cbo  and  define  the  group 
operation  as  composition.  4>fc0  then  forms  a  group,  called  the  holonomy  group  at  6o. 
(Assuming  B  is  connected,  it  is  easy  to  see  that  and  ‘Fftj  are  conjugate  for  any 
two  6o,  6i  G  B.  Thus  if  B  is  connected  we  have  simply  <I>,  the  holonomy  group  of  the 
connection.) 

Definition  2.4.  Given  a  bundle  ir  :  P  — >  B,  a  connection  on  the  bundle,  and  a 
closed  curve  b(t)  in  the  base  space,  the  geometric  phase  is  the  holonomy  along  the 
curve  b(t).  ■ 

The  geometric  phase  can  be  calculated  by  solving  the  ordinary  differential  equa¬ 
tion  given  in  equation  (2.2). 

3.  The  moving  systems  approach  .  Inspired  by  classical  examples  such  as 
the  Foucault  pendulum  and  the  ball  in  a  hoop,  one  is  led  to  consider  the  question 
of  variations  in  phase  space  as  a  parameter  in  a  mechanical  system  is  slowly  varied 
along  a  closed  path  in  parameter  space,  as  in  a  moving  system.  Examples  in  quantum 
physics,  optics,  and  other  settings  [6,  26,  24]  reveal  that  the  essential  calculation  is  a 
geometric  one  and  is  in  fact  captured  by  a  holonomy. 

Marsden,  Montgomery,  and  Ratiu  have  developed  a  modern  geometric  approach 
to  understanding  moving  systems  that  utilizes  the  tools  of  Ehresmann  connections 
on  fiber  bundles  [18].  Here  we  present  a  brief  synopsis  of  their  approach.  Let  S 
be  a  Riemannian  manifold  and  let  M  be  the  space  of  embeddings  of  a  manifold  Q 
into  S.  We  think  of  S  as  the  ambient  space  in  which  Q  is  being  moved  and  of  Q 
as  the  configuration  space  for  a  system  of  interest.  A  tangent  vector  to  M  at  m 
is  a  map  um  :  Q  — >  TS  such  that  um(q)  G  Trn(qyS.  Given  a  tangent  vector  um(q) 
one  can  construct  a  tangent  vector  to  TqQ  as  follows.  Relative  to  the  metric  on  S, 
orthogonally  project  um(q)  to  Tm^q)m{Q)  C  ( Tqm)(TqQ ),  denote  this  vector  u^n{q), 
and  then  pull-back  u^(q)  by  Tm^1  to  TqQ.  This  natural  construction  defines  an 
Ehresmann  connection  on  the  product  bundle  7r  :  Q  x  M  — >  M  as  follows. 

Definition  3.1.  [18]  The  Cartan  connection  on  n  :  Q  x  M  — >  M  is  given  by 
the  vertical-valued  one-form  yc  defined  by 

(3.1)  7 c(q,m)(vq,um)  =  (vq  +  (T~1moull)q,0). 


The  Cartan  connection  induces  a  connection  on  p  :  T*Q  x  M  — >  M  as  follows. 
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Definition  3.2.  [18]  The  induced  Cartan  connection  on  p  :  T*Q  x  M  — >  M 

is  given  by  the  vertical-valued  one-form  "f0  defined  by 

(3.2)  m)  (J-^aq  i  Um)  {Uaq  T  A-p(Um)  (&q)  ?  0) 
where  Vfum)  is  the  function  defined  by 

(3.3)  CP(wm))ag  =  aq  ■  (r'moa^fg) 

and  X-p^Um^  is  the  Hamiltonian  vector  field  ofV(um).  ■ 

To  separate  the  effects  of  the  imposed  motion  on  the  system  (as  defined  by  the 
embeddings  nit)  from  the  nominal  dynamics  (when  the  imposed  motion  is  zero)  we 
use  the  ideas  of  averaging.  Abstractly,  we  assume  we  are  given  a  left  action  of  a  Lie 
group  G  on  T*Q.  The  average  of  the  connection  form  7  is  defined  by 

(3-4)  {!)  =  jG9*{l)d9 

where  dg  is  a  left  Haar  measure  and  |G|  is  the  total  volume  of  G.  From  this  we  have 
the  following  definition. 

Definition  3.3.  [18]  The  Cartan-Hannay-Berry  connection  on  p  :  T*Q  x 

M  — »  M  is  given  by  the  vertical-valued  one-form  7  defined  by 

(3-5)  l{otq,m)(Uaq,Um)  =  ( Uaq  +  A(P(Mm))(a9),0) 

where  (•)  denotes  the  average  with  respect  to  the  action  of  the  Lie  group  G.  ■  In  [18] 
Marsden,  Montgomery,  and  Ratiu  show  that  this  is  an  Ehresmann  connection.  The 
horizontal  lift  of  a  vector  field  Z  on  M  relative  to  7  is 

(3.6)  ( horZ)(aq,m )  =  (-X(p(Z(m)))(aq),  Z(m))  . 

Definition  3.4.  [18]  The  holonomy  of  the  Cartan-Hannay-Berry  connection  is 
called  the  Hannay-Berry  phase  for  a  moving  system.  ■ 

4.  The  rotating,  vibrating  ring  .  In  this  section  we  derive  the  Hannay-Berry 
phase  for  the  vibrating  ring.  Using  the  moving  systems  approach  we  find  an  explicit 
formula  for  the  phase  shift  under  linearizing  assumptions  and  show  that  this  result 
matches  that  of  G.H.  Bryan  [8].  In  this  work  we  are  interested  in  the  effects  of  the 
imposed  rotatory  motion  and  as  a  consequence  choose  to  simplify  the  analysis  of  the 
ring  dynamics  by  assuming  the  ring  has  no  cross-sectional  area.  This  choice  also  allows 
a  direct  comparison  to  the  results  derived  by  Bryan.  A  more  comprehensive  treatment 
based  on  the  geometrically  exact  theory  of  rods  could  be  developed  to  understand  the 
detailed  dynamics  of  the  ring  itself  (see,  for  example,  [20]). 

Consider  a  thin  ring  of  length  L  and  line  density  cr.  The  body  is  given  by  B  = 
{b  :  b  €  [0,  L}}.  Let  6  be  the  mapping  given  by 

9  :  B  ->  S1 


allowing  us  to  parametrize  the  ring  by  9  €  [0,  27t]  .  We  define  the  reference  configu¬ 
ration  to  be  a  circular  ring  of  radius  a  centered  on  an  inertial  reference  frame  and 
( w{9 ),  7 (0))  to  be  the  radial  and  angular  deformations  from  this  reference  respectively. 
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To  maintain  integrity  of  the  ring  we  require  w(0)  =  w(2tt)  and  7(0)  =  7(27r).  In  stan¬ 
dard  cylindrical  coordinates  the  configuration  of  the  ring  is  given  by  (a+w(8),  8+7(8)). 
Let  C  be  the  space  of  all  smooth  deformations  of  the  ring.  (We  do  not  discuss  here 
the  explicit  infinite  dimensional  manifold  structure  for  C  and  associated  structures, 
although  it  is  standard  as  in  [17]).  Since  we  are  interested  in  imposed  rotational 
movements  of  the  ring  (as  a  sensor),  we  split  7 (9)  =  ip  +  a(9)  with  a(0)  =  a(27r) 
where  ip,  independent  of  6,  is  a  global  rotation. 

We  now  use  the  following  argument  of  Rayleigh  [23] .  Since  the  ring  is  thin  the 
forces  resisting  bending  are  small  in  comparison  to  those  which  resist  extension.  In 
the  limiting  case  of  an  infinitely  thin  ring  the  flexural  vibrations  become  independent 
of  any  extension  of  the  circumference  as  a  whole  and  one  may  assume  that  each  part 
of  the  circumference  retains  its  natural  length  throughout  the  motion.  Under  this 
condition  we  say  the  ring  is  inextensible.  Viewing  the  deformed  ring  as  a  curve  in  H2, 
a  point  on  the  curve  is  given  in  Cartesian  coordinates  by 

,  ,  f  x(9)  \_(  (a  +  w(8))  cos(9  +  7(0))  \ 

1  j  ^  »(*)  )~  V  (a  +  w(0))sin(0  +  7(0))  )  ’ 

Equating  the  lengths  of  an  arbitrary  section  of  the  circumference  of  the  reference 
configuration  to  the  length  of  the  same  section  in  the  deformed  configuration  yields 


(4.2) 


Since  8\  and  82  are  arbitrary  we  have 


<9y 


0  -  (a  +  w)  1+wx  +  (  ITT  )  ~a 


88 


dw 


88 


(4.3)  =  2 aw  +  w2  +  2(a2  +  2 aw  +  w2)  +  (a2  +  2 aw  +  w2) 

89 


(9 7 

\88 


From  here  on  we  assume  the  deformations  are  small  and  so  we  keep  only  terms  to 
first  order  in  equation  (4.3).  The  inextensibility  condition  then  requires  that 


(4.4) 


87  da 

W  =  -aM=~aW 


From  the  above  the  space  C  is  given  by 

(4.5)  C  =  {(ip,  a)\ip  €  S'1,  a  :  S1  — >  S1,  a(0)  =  a(27r),  a  smooth}. 

Any  W  £  has  the  form  (with  overdot  denoting  the  partial  derivative  with 

respect  to  time) 


(4.6) 


W  = 


iP 

a(9) 


The  kinetic  energy  is  easily  verified  to  be 


(4.7) 


KE  =  - 
2 


/»27T 


da 
~  88 


■  2ipa(9) 


a2 


da 

~88 


j3dd 


10 


S.  ANDERSSON  AND  P.S.  KRISHNAPRASAD 


where  equation  (4.4)  has  been  used  to  express  w  in  terms  of  a.  This  defines  an  inner 
product  on  C  given  by 


{WUW2)  = 


1  - 


da 


(4.8) 


80 

dd\  da 2 
~d0~d6 


Iip2  +  4>2Ui  (0)  +  ipid2(0)  +  di 
3d0. 


As  in  Bryan,  we  take  the  potential  energy  due  to  the  bending  of  the  ring  to  be 
proportional  to  the  change  in  curvature  squared  of  the  ring.  That  is 

(4.9)  V=^  J  (Ka(0)  -  Ka=0  {0)f  add 

where  na  ( 6 )  is  the  curvature  of  the  ring  at  the  material  point  0  under  the  deformation 
a  and  (3  is  a  material  constant.  For  a  curve  ( r(t),<j>(t ))  defined  in  polar  coordinates, 
the  curvature  is  given  by 


(4.10) 


n(t)  = 


(2  r2(f)  +  rrcj)  —  rrcj)  +  r2cjr) 

(f2  +  r24>2)i 


In  the  0— parametrization,  the  configuration  of  the  ring  under  the  deformation 
(a,  w)  is  given  by  the  curve 


r(0)  =  a  +  w(0) 

4>(0)  =  0  +  ip  +  a(0) 


Using  equation  (4.11)  in  (4.10)  to  express  the  curvature  of  the  ring  under  the  de¬ 
formation  ( w,a ),  we  find  (keeping  only  terms  to  first-order  in  the  numerator  and 
denominator) 


(4.12) 


^a,w(^) 


(— afgr  +  a2  +  2 aw  +  3 a2 1|) 
(a2  +  3a2w  +  3a3|^) 


Writing  the  curvature  in  terms  of  a  alone  using  the  inextensibility  condition  of  equa¬ 
tion  (4.4)  yields 


(4.13) 


Ka(0) 


a 


From  this  we  have  na=o  (0)  =  Using  these  results  in  the  potential  energy  of  equation 
(4.9)  yields 


(4.14) 


/*27T 

d3a  da 

L 

~W  +  ~d0 

dO. 


In  Bryan’s  original  work  the  potential  energy  included  also  a  term  capturing  the 
work  done  in  stretching  the  ring  and  the  work  done  against  an  attracting  force  which 
he  introduced  to  separate  the  effects  of  the  centrifugal  force  from  the  remaining  terms. 
These  two  terms  are  related  by  a  simple  equation  involving  the  rate  of  the  imposed 
rotation  and  at  the  conclusion  of  his  analysis  Bryan  chooses  the  attracting  force  so  as 
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to  cancel  the  tension,  leaving  only  the  work  done  in  bending.  Here  we  take  a  simpler 
approach,  similar  to  Rayleigh,  and  omit  those  terms  at  the  outset. 

The  standard  Lagrangian  function  is  defined  to  be  the  kinetic  minus  potential 
energies  and  is  given  here  by 


L=i 

2 


da 
~  d9 


2  ipa{9)  +  of 


(4.15) 


A 

2  a  . 


d3a  da 

d(P  +  m 


i  2 


dd. 


Consider  now  the  following  action  of  S1  on  C. 


(4.16) 


$g(ip,a)  =  (ip  +  g,a). 


We  have  the  following  lemma. 

Lemma  4.1.  (C,  (•,•),  H,  S'1)  is  a  simple  mechanical  system  with  symmetry  where 
the  action  of  S1  on  C  is  given  by  equation  (4.16).  (For  a  definition  and  discussion  of 
simple  mechanical  systems  with  symmetry  see  Section  4.5  of  [1].,) 

Proof.  Immediate  since  both  (•,•)  and  V  are  invariant  under  the  given  action  of 
S1  on  C.  □ 

Since  the  given  action  is  both  free  and  proper,  the  reduced  space  Q  =  C/S1,  given 
by 

(4.17)  Q  =  {a  :  S1  — +  S1|a(0)  =  a(27r),  a  smooth} 


is  also  a  manifold.  To  fix  notation  in  relation  to  the  general  theory  presented  in 
Section  3,  we  note  that  Q  =  C/S1  is  the  configuration  space  for  the  ring  and  S  =  C  is 
the  ambient  space  in  which  Q  is  moved.  To  slowly  rotate  the  ring  we  set  ip  =  ip0  +  fit 
(identifying  ip  =  0  with  ip  =  27r)  for  some  small  Cl  and  some  fixed  initial  offset  ipo  so 
that  the  embedding  from  Q  to  S  is  given  by 

(4.18)  mt{a{9 ))  =  ( ip0  +  fit,  a(9)). 


4.1.  The  nominal  dynamics.  The  nominal  dynamics  is  given  by  setting  fl  =  0 
in  equation  (4.18).  Applying  this  to  equation  (4.15)  yields  the  nominal  Lagrangian. 


L 


Nom 


(«,*)=/  ^ 


i  da  \  • 


•2 


2 

VM 


(4.19) 


A 

d3a  da 

2\ 

2  a 

J93  +  09  _ 

1 

d9 


The  action  integral  for  this  Lagrangian  is  defined  to  be 

rb 


(4.20) 


J(a,  a)  =  f  LNom(a,a)dt. 
J  a 


The  Euler-Lagrange  equations  for  this  system  are  found  by  applying  Hamilton’s 
principle  of  critical  action  (see,  e.g.  [1])  which  states  that 


(4.21) 


r-b 

6J(a,  a)  =  6  LNom{a,  a)dt  =  0 
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for  all  variations  among  paths  77(f)  in  Q  with  fixed  end-points.  Applying  the  variation 
yields 


8J(a,  a)  =  — 

d 

de 


J(a  +  erj,a  +  erf) 


€=0 

pb  p2ir 

g— 0  J  &  JO 


era 

2 

2  a 


1  - 


9  (a  +  erf) 

00 


(a  +  e?))2  + 


9(d  +  erf) 

86 


d3(a  +  e?7)  |  8{a  +  erf) 

W3  +  86 


dddt 


nb  />2lT 

la  Jo 


da  \  .  2  dij 


86 )  86 


da 


—  (  1  ^7T  )  ^  ^71  +  (  1  T 


da  8i 7 


/? 


93a  9a 

863  +  86 


83ri  8i  7 

903  +  90 


90 

d0cif. 


96*  86 


Using  integration  by  parts  repeatably  on  the  space  variable  and  the  fact  that  for  any 
element  £(•)  £  Q  we  have  £(0)  =  £(27t),  the  variation  of  the  action  can  be  rewritten 
as 


b  r2ir 


8J(a,  a)  = 


a  JO 


8a\  da 


82c 


2  1  -  1177  -Zi “  v  +  1  -  ivy  ar]  -  yryyf? 


1 

a 


86  /  86 


86a  94a  ®2<X 

906  +  2  ~d¥  +  W 


862 


77  >  dddt. 


da  V 


82c 


9<2 ; 


96»2 


Using  integration  by  parts  once  again,  this  time  on  the  time  variable,  and  utilizing 
the  end  point  condition  on  the  variations  77  yields 


b  (‘2-k 


SJ(a,  a)  = 


a  JO 


4  1- 


8a\  da  .  82a 

902 


86  /  86 


86a  94a  ®2a 

906  +29(U  +  -QQ2 


a2-  1- 


8a 

90 


92d 

902 


r)d6dt. 


Since  this  must  equal  zero  for  all  variations  77  we  obtain  the  Euler-Lagrange  equations 
for  the  nominal  system. 


(4.22) 


0 


era 


92  a 
902 


a 


f3  86a  94a  ^2<x 

+  a  906  +2904  +  902 


To  simplify  this  difficult  nonlinear  partial  differential  equation  we  use  the  assump¬ 
tion  that  the  deformations  are  small  and  replace  the  above  equation  by  its  lineariza¬ 
tion,  resulting  in  the  following  equation  of  motion  for  a. 


(4.23) 


82a 

" 

(3 

’  96  a 

84a 

82a 

862 

—  a 

+  - 
a 

892 

+  2W  + 

862_ 

At  this  stage  we  state  our  intention  to  do  all  of  the  calculations  associated  to 
the  application  of  the  moving  systems  approach  to  the  vibrating  ring  problem  in  a 
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convenient  set  of  coordinates,  namely  the  coeffecients  of  a,  d  expressed  in  a  Fourier 
basis.  We  first  express  the  nominal  dynamics  in  these  coordinates  and  in  the  following 
section  do  the  holonomy  calculations  in  the  same  coordinates  (after  truncation  of  the 
Fourier  series).  In  these  coordinates  a  has  the  form 

OO 

(4.24)  a(9)  =  ^  [ Ak  cos (kd)  +  Bk  sin(/d?)]  . 

fc= 1 

The  deformation  a{9)  is  not  allowed  to  contain  any  global  rotations  and  so  the  con¬ 
stant  coefficient  is  set  to  0.  Inserting  this  expression  for  a  into  the  equation  of  motion 
(4.23)  results  in  the  equation 


0  =  ^  |cra3(l  +  k2)  Ak  cos {k6)  +  Bk  sin (kO) 

fc= l 


(4.25) 

+  — (fc6  —  2 fc4  +  k2)  [Akcos(kd)  +  Bk  sin(fc<?)]  |  . 

Collecting  terms  in  cos(0)  and  sin(0)  and  setting  them  separately  to  zero  gives  the 
following  set  of  ordinary  differential  equations  for  the  Fourier  coefficients. 

(4.26) 

a  &  k2(k2  -  l)2  a  2  a 

Ak~  aai  k2  +  i  Ak~  pkAk, 

(4.27) 

p  P  k2{k2  —  l)2  a  2 

Bk~  aa4  k2  +  1  Ak~  9kBk 

which  defines  the  frequencies  rjk-  This  result  is  in  agreement  with  a  derivation  of 
Rayleigh  [23]  and  defines  for  each  k  a  pair  of  uncoupled  oscillators  with  common 

frequency  rjk- 

The  solution  to  this  system  is  given  by 

(4.28) 

_  Ak 

Ak(t)  =  Akcos(r]kt)  H - sm(r]kt), 

m 

(4.29) 

^  B 

Bk(t)  =  Bk  cos (rjkt)  H - sin (r]kt) 

Tli 

Vk 


where  Ak,  Ak,  Bk,  and  Bk  are  given  by  initial  conditions.  The  Hannay-Berry  phase 
is  defined  as  the  holonomy  on  a  trivial  bundle  involving  the  cotangent  space  of  the 
system.  It  will  prove  useful,  then,  to  have  the  time  evolution  of  the  conjugate  momenta 
for  the  nominal  system.  By  inserting  the  Fourier  expansion  for  a  into  the  nominal 
Lagrangian,  equation  (4.15),  and  applying  the  Legendre  transform  we  obtain 

dL 

PAk  =  — y  =  (1  +  k2)aa3TrAk, 
oAk 

dL  o  o 

PBk  =  =  (!  +  k2)aa3TrBk. 

oBk 

Thus  the  solution  to  the  nominal  system  expressed  on  the  cotangent  bundle  is  given 

by 

(4.32)  Ak(t)  =  Ak  cos(rjkt)  +  P^k  - sin (r]kt), 

(1  +  k^jaa^mjk 


(4.30) 

(4.31) 
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PAk(t)  =  -(1  +  k2)aa3npkAk  sin (r)kt)  +  pAk  cos(r]kt, ) 


Bk(t )  =  Bkcos(pkt) 


(1  +  k2)aa3nrik 


sin(%f), 


(4.35)  PBk(t )  =  -(1  +  k2)aa3npkBk  sm(pkt)  +  PBk  cos(r]kt). 

4.2.  The  Hannay-Berry  phase.  From  equation  (4.18),  the  velocity  vector  of 
the  motion  in  S  is 

(4.36)  4  (mt(a(9)))  =  (0,  d(9))  +  (Cl,  0). 

From  this  we  see  that  the  vector  field  corresponding  to  the  imposed  motion  is 

(4.37)  Z(mt(a(9)))  =  (  q  )  ' 

For  ease  of  notation  define  Z  =  Zt(mt(a(9))) .  The  projection  of  this  tangent  vector 
to  Tmt^mt{Q)  with  respect  to  the  metric  of  S  is  given  by  ZT  =  Z  —  Z1-  where 
(Z-1,  A)  =  0  MX  G  Tmt(qpTnt{Q) .  From  equation  (4.6),  any  vector  X  G  Tmt^mt{Q) 
has  the  form 


where  Y  G  TqQ.  Then,  from  equation  (4.37)  and  the  fact  that  ZT  G  Tmt(9)?nt((5),  we 
can  write  ZT  =  Z  —  Z1-  as 


for  some  Yzt,Z±2  G  Tq(Q)  and  Z3-1  G  T^S1.  From  equation  (4.39),  Z3-1  =  O. 
Applying  equation  (4.8),  the  orthogonality  condition  states 

0=((F!,ZJ-2),(0,F)) 


('"iO  ^Y  +  zUY1  +  m^e 


In  what  follows,  we  express  the  orthogonality  condition  of  equation  (4.40)  in  the 
Fourier  basis  and  thus  derive  an  explicit  formula  for  Z ±2  (see  (4.51)  and  (4.52)  below). 
Using  the  Fourier  series  representation  for  a,  a  tangent  vector  Y  G  TqQ  has  the  form 

OO 

(4.41)  Y  =  ^2  cos (kO)  +  YBk  sin(A;0)] 

k= 1 

Adopting  the  notation  c (9)  =  cos (9)  and  s(9)  =  sin(0),  the  orthogonality  condition  is 
given  by 

„2tt  /  oo  \  2  /  oo 

0=  /  1  -  k  [Bkc(k9)  -  Afcs(fc0)]  Akc(k9)  +YBks(k9)] 

\  k — 1.  J  \  k — 1 

oo  oo  \ 

+  E  E  [ZAAkd)  +  ZbA™)]  1YaA19)  +  YbAW)} 

fc= 1  i= 1  / 

oo  oo 

(4.42)  +  E  E  kl  (ZbA™)  -  Z%*m)  (YbAW)  -  YaAW))  ™*d9. 
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Using  the  following  identities 


c(k9)c(l9)d9  =  /  s(kd)s(W)de  =  nSki, 


c(k9)c(l0)c(m0)c(n6)d9  =  —Skimn, 

K  t 

s(k9)s(l9)s(m9)s(n9)d9  =  — Skimn , 

r  7 r 

c(k9)c(l9)s(m9)s(n9)d9  =  —  Skimn 


where 


(4.47)  5kl  = 


1  if  k  =  1 
0  otherwise 


dklmn  — 


1  if  k  =  l  =  m  =  n 
0  otherwise 


and  the  fact  that  all  other  combinations  of  sin  and  cos  appearing  in  equation  (4.42) 
integrate  to  0  over  [0,  2ir\,  the  orthogonality  condition  can  be  reduced  to 


0  =  [ (z\l  +  ^A\Zil  +  B 2kZ%  -  ^ AkBkZi 


-2 ClkBk  +  k2z\l}  YAk 


+  iJTAlzBl  +  T bIzbI  - 


+2 nkAk  +  k2Z^YBk  . 


This  holds  for  every  Y  €  TqQ  and  so  for  every  k  we  have 


(4.49)  Z^l  +  -A2Z^k  +  —B2Z\l  -  -AkBkZ%  -  2 ClkBk  +  k2Z £  =  0 


(4.50)  Z 


**  4 

k2 


I  K  ^-k  2  nk 

jBlZil  -  kdA„B„Zil 


r2  +  —  Alzi2  +  —Bizi2  -  —AkBkZi2  +  2 ClkAk  +  k2Zi2  =  0. 


Solving  these  coupled  equations  for  Z\2  and  Z^2  yields 


Zi2  =  -Cl 


2fc(l  +  k 2  +  ^fA2k  +  ^Bg)Bfc  -  k3A2Bk 
Dk(Ak,  Bk ) 

~2fc(l  +  fc2  +  +  ^Bg)Afc  -  k3AkBl 

Dk(Ak,  Bk) 


where 


Dk(Ak,  Bk)  —  (1  +  k2  +  —A I  +  -^-Bk)(  1  +  k2  +  — j--4fc  +  ~^Bk) 
~  —  A2B2 

^  Aknk 


(4.53) 
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and  so,  representing  the  tangent  vector  by  its  coefficients  at  each  k, 

(  \ 


(4.54) 


z1-  = 


V 


O 


-tt 


2k(l+k2  +  ^Al  +  ^rBl)Bk-k3AlBk 
Dk(Ak,Bk) 


k= 1 


Dk(Ak,Bk) 


k=l  / 


Inserting  equation  (4.54)  into  equation  (4.39)  gives  the  projection  of  the  tangent  vector 
of  the  imposed  motion  onto  Tmt^mt(Q)  to  be 


/ 


(4.55) 


= 


V 


-fi 


n 


o 

+. 1 

_ '  4 _ 

Dk(Ak,Bk ) 

2fc(l +k2+’£A2k+*£-Bl)Ak-k3AkBl 
Dk(Ak,Bk ) 


\ 


j  j  fc=i 


C  =  1  ) 


The  pull-back  of  Z1  to  TqQ  by  [Tm\  1  is  given  by 


A 


(4.56)  Z(q)  =  [Tm]~1Z1  = 


l 


~n 


o 


2k(l+k2+^-Aj  +  SrB2k)Bk-k-sA2kBk 
Dk(Ak,Bk ) 

2k(l+k2  +  ^Al  +  ^Bl)Ak-k3AkBl 
Dk(Ak,Bk) 


k= 1 


fc=l 


) 


where  Z(g)  is  defined  for  ease  of  notation.  Recalling  that  the  deformations  are 
assumed  to  be  small,  the  above  expression  is  expanded  in  a  Taylor  series  about 
Ak  =  Bk  =  0  V  k  and  only  the  first  order  terms  kept.  This  yields 


(4.57) 


Z(q)  = 


/  {^BC, ) 


l  {SMC.  / 


Let  qk  =  (Ak,Bk)  so  that  the  coordinates  on  Q  are  q  =  {qkjkLi-  The  conjugate 
momenta  can  then  be  written  pk  =  ( PAk,PBk )•  Define 


(4.58) 


ry  /  \  A  2  k  ft 

Zk[qk)  = 


0  -l 

l  +  k2  V  1  0 


Ak 

Bk 


2  kCl 
l  +  k2 


S(l)qk 


which  also  defines  the  skew  symmetric  matrix  5(1) .  With  this  definition  the  projected 
vector  can  be  expressed  as  Z(q)  =  {Zk(qk)}kLi-  To  avoid  technical  difficulties  we 
assume  that  only  N  Fourier  modes  are  active  where  N  is  some  positive,  finite  number. 
Using  this  assumption,  the  function  V(Z)  defining  the  horizontal  lift  relative  to  the 
induced  Cartan  connection  (as  in  equation  (3.3))  is  given  by 


N  r 


(4.59) 


V(Z)(q,p)=p.Z(q)  =  Y/ 


fc= l 


2  kn 


1  +  k2 


Pk  ■  S(l)qk. 
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Define  Qk  =  {( Ak,Bk )  €  IR2}  so  that  Q  =  U k=1Qk-  The  configuration  space  is 
then  the  Cartesian  product  of  N  copies  of  R2  and  each  coordinate  qk  and  conjugate 
momenta  pk  can  be  identified  with  a  vector  in  R2.  Extend  these  vectors  to  R3  by 
letting  the  third  coordinate  of  each  be  zero.  Let  x%k  be  a  unit  vector  at  the  origin  of 
the  kth  copy  of  R2  along  this  third  direction.  With  these  identifications,  in  equation 

(4.59)  we  replace  pk  ■  S(l)qk  by  (qk  x  pk)  ■  x3k.  Define 

(4.60)  4  =  (qk  x  pk)  ■  x3k 
so  that 

JL  2  k 

(4-61)  V(z)(q,p)=^2—¥mk. 

k=l  1  ' 

Let  T  be  the  subset  of  C°°  functions  on  T*Q  defined  by 

(4.62)  T  =  ^f(q,p)\f(q,p)  =  ^  akfk(Ak,  Bk,pAk  ,pBk),  ak  G  TR,  fksmooth^  . 

Since  the  time  solution  to  each  Fourier  mode  for  the  nominal  system  is  periodic,  we 
can  define  an  average  on  T  with  respect  to  the  flow  by 

N 

(4-63)  (/)  =  '52ak(fk)k 

k=  1 


where 

(4-64)  (fk)k  =  ^  jjk  ($)*  fkdt 

where  is  the  flow  corresponding  to  the  kth  Fourier  mode  of  the  nominal  system. 
From  equations  (4.60)  and  (4.61)  it  is  clear  that  P(Z)  €  T .  We  have  the  following 
useful  lemma. 

Lemma  4.2.  Ik  is  constant  along  the  trajectories  of  the  nominal  system. 

Proof. 

dlk  d  .  \  ^ 

^-  =  S(®X»)-X  3, 

=  Jt(pBkAk-pAkBk) 

=  i>Bk  Ak  +  pBk  Ak  -  pAk  Bk  -  pAk  Bk 
=  (1  +  k2)cra3Trqk  ^ BkAk  +  BkAk  —  AkBk  —  AkBk ''j 

=  (1  +  k2)cra3nr]k  ( BkAk  -  AkBk ) 

where  in  the  second  to  last  step  we  have  used  the  definition  of  the  conjugate  momenta 
for  the  nominal  system  in  equations  (4.30,4.31).  From  equations  (4.26,4.27),  along 
the  trajectories  of  the  nominal  system  we  have 

( BkAk  —  AkBkJ  =  —r/2  ( AkBk  —  AkBk)  =  0. 


(4.65) 
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□ 

The  average  of  V(Z)  over  the  nominal  dynamics  is  then 

*  2k  i,  2  k 

(4.66)  (V(Z))( ,.p)  =  •£  =  X 

k— 1  fc— 1 

Noticing  that  this  function  depends  only  the  coordinates  (■ qk,Pk )  through  Ik,  we  move 
to  the  coordinates  for  the  averaged  dynamics  in  phase  space  defined  by  (fa,  Ik,  Pk,Ppk ) 
where  pk  =  ((A|  +  B2))^ .  Here  fa  is  conjugate  to  Ik  and  pPk  is  conjugate  to  pk- 
With  these  coordinates  the  horizontal  lift  of  Cl  relative  to  the  Cartan-Hannay- 
Berry  connection  (as  in  equation  (3.6))  is  given  by 

(«7)  (-A,, *<„»,«)  =  ({-rfp^L  ■«)  ■ 

The  geometric  phase  is  neatly  split  into  a  phase  change  in  each  fa  independently.  If 
the  loop  in  M  is  parametrized  by  t  £  [0,  T]  where  T  =  is  the  time  to  complete  one 
full  revolution  of  the  ring,  then 

r T  2k  2k 

(4.68)  A  fa  =  ~Jo  YT¥ndt  =  ~l  +  ¥nT=~2lT 

After  one  full  rotation  of  the  ring  each  vector  qk  has  been  rotated  by  the  angle  A  fa. 
In  practice  one  expects  to  use  only  one  mode  under  a  resonant  drive  (as  in  [21,  22]). 

Note  that  we  perform  an  average  with  respect  to  the  flow  of  the  nominal  dynamics 
on  a  special  class  of  functions  as  in  equation  (4.62).  This  agrees  with  a  group  (here 
S1)  average  as  in  Section  3  (Definition  3.3)  when  we  restrict  the  nominal  dynamics 
to  a  single  mode. 

5.  A  comparison  with  the  results  of  Bryan  .  In  [8]  Bryan  uses  classical 
variational  techniques  to  derive  the  equations  of  motion  for  a  thin  ring  of  radius  a 
undergoing  a  steady  rotation  about  its  central  axis  with  angular  velocity  Cl.  His 
analysis  uses  two  polar  coordinate  systems,  the  first  fixed  in  space  and  the  second 
rotating  with  angular  velocity  Cl.  If  in  the  undeformed  state  the  coordinate  systems 
are  given  by  (a,  </>)  and  (a,  9)  we  have 

(5.1)  <j>=e  +  cit 


2k 

1  +  fc2 


and  9  is  constant  for  any  particle  of  the  ring.  Let  the  tangential  and  radial  displace¬ 
ments  of  a  particle  of  the  ring  be  given  by  v  and  w  respectively  so  that  the  new 
polar  coordinates  are  (a  +  w,  <fi  +  v/a)  and  (a  +  w,  9  +  v/a )  in  the  two  systems.  The 
deformations  v,  w  are  assumed  to  be  small.  The  assumption  of  inextensibility  yields 


(5.2) 


dv 


w  =  — 


d9 


as  before.  As  discussed  in  the  previous  section,  Bryan  includes  work  done  against 
the  tension,  T,  to  stretch  the  ring  and  against  an  attractive  force,  p.  To  match 
his  derivation  with  the  model  we  have  chosen,  we  set  these  terms  to  zero.  Taking 
variations  on  the  total  energy  and  setting  them  to  zero,  we  find  the  following  equation 
of  motion 


..  d2v  dv 

~  V  89 2  ^  89 


Clz 


,82v 

W 


(3  82 

aa4  892 


(5.3) 


82 

892 


2 


V. 
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Bryan  then  assumes  that  the  deformations  are  of  the  form 
(5.4)  v  =  cos(kd  +  pt). 


Inserting  this  into  the  equation  of  motion  yields  the  following  two  solutions. 


(5.5) 

(5.6) 
where 


v  =  Acos  (  kd  + 


2k 


1  +  k2 


v  =  Acos  kd  + 


2k 


1  +  k2 


O t  T 


fit  —  U>kt 


(5.7) 


2  _  p  fc2(fc2~  l)2  _  fl2fc2(fc2  -  3) 
k  era4  1  +  k2  (1  +  k2)2 

2  fl  2k2(k2  —  3) 

=  Vk  (1  +  k2)2  ' 


Notice  that  by  retaining  the  terms  in  fi2  there  is  a  slight  decrease  in  the  frequency 
of  vibration  from  r^..  This  can  be  understood  as  a  “softening”  of  the  material  and 
corresponds  to  the  spurious  softening  that  occurs  in  the  theory  of  rotating  rods  if 
the  models  are  linearized  prematurely,  that  is  before  the  effects  of  external  rotation 
are  considered.  The  geometrically  exact  theory  handles  this  issue  properly  and  we 
believe  that  extending  the  model  of  the  ring  using  this  theory  will  prove  useful.  (We 
note  also  that  if  the  attractive  force  p  introduced  by  Bryan  is  kept  in  the  equations 
of  motion  and  set  equal  to  O2  then  the  resulting  system  shows  an  increase  in  the 
frequency  of  vibration  due  to  the  imposed  rotation.  See  [13]  for  comments  on  similar 
ad  hoc  methods  in  the  theory  of  rotating  beams.) 

Assuming  the  amplitude  A  in  the  two  solutions  is  the  same,  the  final  solution  in 
the  fixed  frame  is  given  by 


(5.8) 


v  =  2  Acos(ukt)cos 


k2  —  1 
k2  +  1 


Bryan  then  recognizes  that  this  corresponds  to  an  oscillation  with  2k  nodes  where  the 
position  of  the  nodes  precess  in  retrograde  around  the  ring  with  angular  velocity 


(5.9) 


k2  —  1 
k2  +  1 


n. 


If  we  write  this  in  the  rotating  system  we  have 

<5-10>  (Srr  -  0  n  =  -FTTn 


and  after  one  rotation  the  nodes  have  precessed  in  the  moving  frame  by 
(5-u) 


To  compare  this  to  our  results  in  Section  4  we  must  first  restrict  our  solution  to 
a  single  mode  so  that 


(5.12) 


a(9,t )  =  Ak{t)  cos  (kO)  +  Bk{t)  sin  (kd). 
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When  ( Ak,Bk )  is  viewed  as  a  vector  in  IR2,  the  effect  of  the  geometric  phase  is  seen 
to  be  a  rotation  of  this  vector  about  the  origin  where  the  counter-clockwise  direction 
is  defined  to  be  positive.  Using  equation  (4.68),  the  rotated  vector  at  the  end  of  one 
revolution  of  the  ring  is  given  by 


Ak(T) 

Bk{T) 


sin 


(5.13) 


~27T^)  -sin(-27r4|i) 

cos  (”27rw) 

Afc(0)  COS  -  Bk  (0)  Sill  (~2lT 

Afc(0)  sin  (-27Tj||2 )  +  Bk{ 0)  cos 


Ak{  0) 
Bk{  0) 


Inserting  this  into  equation  (5.12)  and  simplifying  we  get 

2 


(5.14) 


a(6,  T)  =  ^U(O)  cos 

+Ufc(0)  sin  (k 


0  +  2ir 


9  +  2ir 


1  +  k2 

2 


1  +  k2 


which  is  course  the  same  solution  with  the  nodes  rotated  by  —  27t 
with  Bryan. 


l+k2 


agreeing 


6.  Nonlinear  corrections  .  We  now  turn  to  an  investigation  of  corrections  to 
the  geometric  phase  based  on  the  nonlinear  terms  in  the  vector  field  Z{q)  given  in 
equation  (4.56).  It  is  worth  noting  that  these  arise  due  to  the  configuration-dependent 
quadratic  form  defining  the  kinetic  energy.  We  proceed  by  keeping  higher-order  terms 
in  the  Taylor  expansion  of  Z{q).  The  second-order  terms  in  this  expansion  can  be 
shown  to  be  zero.  To  third  order  the  vector  field  is 


(6.1) 


z(q)  = 


/ 

2  kBk 

k2(2+k)A2kBk 

9  k3  Bk  ' 

Y  \ 

1+k 2 

(1+k2)2 

(1+k2)2 

J  k=l 

V 

jo 

2  kAk 

k2(2+k)AkB2k 

9k3  A3  1  1 

N 

1+k2 

(1+k2)2 

(l+fc2)2J  J 

fc=i  / 

Define  the  matrix  U{qk)  by 

(6.2)  U{qk)  = 
With  this,  equation  (6.1)  can  be  written  as 

Z{q)  = 

(6.3) 


Al 

0 


Bl 


2k.n  k2{2  +  k)fl 

S(l)qh  -  ,  ,9x9  U(qk)S{l)qk 


1  +  k2 


{1  +  k2)2 

9  k3f2 


N 


2S(l)U{qk)qk  \ 


(1  +  k2) 

and  the  function  defining  the  Hannay-Berry  phase  is  given  by 


)  fc=i 


N  r 


V(Z){q,p)  =  J2 


2kn  cn\  k2{2  +  k)n 

iPk  ■  b(l )qk  ~  n  ,  ,2,2  Pk  ■  U{qk)b{l)qk 


fc=i  L 


1  +  k2‘ 


{1  +  k2)2 
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(6.4) 


=  £ 


2  kil 

TTfc24 


9  k3n 
(1  +  k 2)2 


■Pk  ■  S(l)U(qk)qk 


k2(2  +  k)il  TT/  ^ 

(1  +  fc2)2  Pk  ■  U(qk)S(l)qk 


9  k3n 
(1  +  k2)2 


pk  ■  S(l)U(qk)qk 


where  in  the  second  step  the  definition  of  Ik  from  equation  (4.60)  has  been  used. 
To  determine  the  Hannay-Berry  phase  we  need  to  find  the  average  of  equation  (6.4) 
over  the  nominal  dynamics.  In  Section  4  we  have  shown  that  the  first  term  in  square 
brackets  in  the  sum  in  equation  (6.4)  is  a  constant  along  trajectories  of  the  nominal 
system.  The  second  and  third  terms,  however,  are  not  constant  and  their  averages 
need  to  be  explicitly  calculated.  For  the  second  term  the  average  of  the  kth  element 
of  the  sum  is  given  by 

,32L 

(Pk  ■  U(qk)S(l)qk)k  =  ^  J  k  [—pAk  (t)  Ak (t) Bk (t) 

(6.5)  +pBk(t)Ak(t)Bl(t)]  dt. 

Using  the  solution  to  the  nominal  system  given  in  equations  (4.32-4.35)  and  making 
a  change  of  variables  in  the  integration,  this  becomes 


1  f2n 

( Pk  ■  U(qk)S(l)qk)k  =  —  J  {  [(1  +  k2)aa?ivqkAk  sin (t)  -  pAk  cos(t)] 


Ak  cos (t) 


pAk  sin(t) 


Bk  cos (t) 


(6.6) 


(1  +  k2)aa37rr]k  J 

+  [-(1  +  k2)aa3TrqkBk  sin(t)  +  pBk  cos(t)] 
pAk  sin(t) 


Ak  cos (t) 


(1  +  k2)aa3nrik 


Bk  cos  (t) 


pBk  sin (t) 

(1  +  k2)aa3irrik 

PBk  sin  (t) 


(1  +  k2)aa3nr]k 


dt 


where,  through  a  standard  abuse  of  notation  in  averaging,  the  hats  have  been  dropped 
on  the  initial  conditions.  The  integration  identities  in  equations  (4.43  -  4.46)  can  be 
used  to  write  the  above  expression  as 

(pk  •  U(qk)S(l)qk)k  =  i  [(AlpBk  -  B\pAk  -  A2kBkpAk  +  AkB2pBk) 

Akp2AkpBk  -  BkpAkp2Bk  -  Bkp3Ak  +  Akp3Bk 
((1  +  k2)aa3Trr]k)2 

Now 


AkPBk  —  BkpAk  —  AkBkpAk  +  AkBkpBk  —  ( AkpBk  —  BkpAk)(Ak  +  Bk) 
(6.8)  =Ik(A2+B2) 

and 


AkPAkPBk  -  BkpAkp2Bk  -  Bkp\k  +  Akp%k  =  ( AkpBk  -  BkpAk){p2Ak  +p%k) 

=  h{p2Ak+P2Bk) 


(6.9) 


22 


S.  ANDERSSON  AND  P.S.  KRISHNAPRASAD 


and  so  equation  (6.7)  takes  the  form 
(6.10)  (pk  •  U{qk)S{l)qk)k  —  vr  I  Aj,  +  Bk  + 


4  l  A  2  ,  D2  ,  Pa,  +p|. 


((1  +  k2)aa3Tnjk)2  )  ' 


Following  the  same  procedure,  the  average  of  the  kth  element  in  the  sum  of  the 
third  term  is 


(Pk  ■  S{l)U{qk)qk)k  =  ^  J  ( PBk(t)A3k{t )  -  pAk{i:)B3k{t))dt 

=  g  [{AkBkpBk  -  A2kBkpAk  +  A\pBk  -  B\pAk ) 

,  -  JSfepl  +  Akp2A  pBk  -  BkPAkp2Bk 


p2ir 


(6.11) 


^  +  PA*+PB« 


((1  +  k2)aa3irpk)2 

,2  i  „2  \ 


o  y  fc  k  '  ((1  +  k2)aa37rVk)2  J  • 
Using  equations  (6.10,6.11),  the  average  of  P(Z)  is  given  by 

N 

(V(Z)((q,p)  = 


fc=i 


2fc 


1  +  fc2 


(6.12) 


/e2  +  14/c3  (  Al  +  Bl+  ((1+fc2fc)o.a3^fc)2 

2(1  +  /c2)2 


From  equations  (4.32,4.34)  we  see  that  the  term  in  parentheses  is  the  average 
of  (A2  +  Bk)  over  the  nominal  dynamics.  We  move  to  the  averaged  coordinates 
(<f>k,  Ik,  pk,pPk)  as  in  the  comments  following  equation  (4.66).  In  these  coordinates, 
the  average  of  V(Z)  has  the  form 


(6.13) 


JV 

(V(Z)((q:p)  =  J2n  h 

fc= l 


2  k 

1  +  k2 


k2  +  14fc3  2 
2(1  +  k2)2  Pk 


The  lift  to  third-order  of  Q  with  respect  to  the  Cartan-Hannay-Berry  connection 
is  given  by 


(6.14) 


2k 

1  +  k2 


f^dkpk 


Uk3 


2(1  +  k2) 


2 

2  Pk 


d 

d(j>k 


k 2  +  14  k3 

5  1 

_  (1  +  fc2)2  _ 

dpPk  j 

N 


fc= 1 


From  equation  (6.14)  we  see  that  both  Ik  and  pk  are  constant.  Thus 


(6.15) 


A  <t>k 


0  T  2 k  k2  +  14fc3 
[l  +  fc2  “  2(1  +  k2)2 
2k  k2  +  14A:3  2' 
ITfc2  ~  2(1  +  k2)2Pk_ 


Pk 


dt 
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Notice  that  the  third-order  terms  act  to  reduce  the  rate  of  nodal  rotation  and  thus 
the  sensitivity  of  a  vibrating  ring  gyroscope  cannot  be  increased  by  increasing  the 
amplitude  of  vibration  and  using  the  nonlinear  effects. 

In  contrast  to  the  earlier  calculation  where  we  kept  in  Z  only  the  terms  linear  in 
configuration  variables,  in  the  present  nonlinear  setting  the  imposed  rotation  causes 
not  only  a  precession  of  the  nodes  of  vibration  but  also  a  drift  in  the  momentum 
conjugate  to  pk  ■  In  practical  devices  the  ring  is  driven  into  a  single  mode  of  oscillation 
and  the  imposed  rotation  sensed  by  measuring  the  drift  rate  of  the  nodal  points  of  the 
vibration.  Thus  the  effect  of  the  second  term  in  equation  (6.14)  will  be  compensated 
for  by  the  drive  electronics. 

It  is  interesting  to  ask  how  large  the  nonlinear  effect  on  the  drift  rate  of  the  nodal 
points  of  the  vibrations  is  in  a  typical  device.  The  micromachined  ring  gyroscope 
of  Putty  and  Najafi  [22]  utilizes  a  ring  of  radius  a  =  500  pm  placed  into  elliptical 
vibration  so  that  k  =  2  with  a  radial  deformation  amplitude  of  0.15  pm..  From 
equations  (4.4,4.24),  the  radial  deformation  for  this  ring  is 

(6.16)  w{9)  =  2 a  [ A2  sin(20)  —  B2  cos(20)] . 


Let  t  =  0  be  the  time  at  which  the  maximum  radial  deformation  is  attained  and 
let  8  =  0  be  the  location  on  the  ring  of  the  maximum  radial  deformation  at  time 
t  =  0.  From  these  definitions  we  have  the  initial  conditions  A2  =  pa2  =  Pb2  =  0  and 
B 2  =  —  2 (so (pirn)  •  Inserting  these  values  into  equation  (6.15)  yields 


(6.17) 


A(f>2  =  —2n 


4 

5 


-  2.61  x  1CT8 


For  the  normal  operation  of  this  device,  then,  the  nonlinear  effects  are  seven  orders 
of  magnitude  smaller  than  the  first-order  terms. 

7.  Conclusions  and  future  work  .  In  this  paper  we  have  applied  the  moving 
systems  approach  developed  by  Marsden,  Montgomery,  and  Ratiu  to  the  rotating, 
vibrating  ring  and  showed  that  the  precession  of  the  nodes  with  respect  to  the  ring  can 
be  understood  as  the  Hannay-Berry  phase.  We  have  also  showed  that  under  the  same 
linearizing  assumptions,  our  result  matches  that  of  Bryan’s  original  analysis.  Using 
the  inherently  nonlinear  nature  of  the  moving  systems  approach  we  then  calculated 
the  effect  of  the  imposed  rotation  on  the  nodal  precession  to  third-order.  These  higher- 
order  terms  serve  to  reduce  the  sensitivity  of  the  device  and  we  therefore  conclude 
that  the  best  performance  of  these  devices  is  obtained  by  operating  them  in  the  linear 
regime. 

There  is  an  underlying  assumption  that  the  rate  of  rotation  of  the  ring  is  slow 
with  respect  to  the  frequency  of  vibration.  Under  this  assumption  the  character  of  the 
nominal  dynamics  is  unaffected  by  the  imposed  motion  (in  the  linearized  setting)  and 
the  effect  of  the  rotation  is  completely  captured  by  the  nodal  precession.  At  rotation 
rates  that  are  not  small  with  respect  to  the  rate  of  vibration  the  dynamics  are  severely 
affected  and  these  techniques  do  not  apply.  However  even  at  slow  rotation  rates  the 
effect  on  the  vibration  is  non-zero.  In  a  future  paper  we  will  report  on  recent  work  in 
which  we  extend  the  theory  of  moving  systems  using  Hamiltonian  normal  form  theory 
to  account  for  the  non-adiabatic  nature  of  the  imposed  motion  (see  also  [2]). 
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